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Wave  spectral  parameter 
Area 

Wave  spectral  parameter 
Composite  (subscript) 

Slamming  parameter 
Froude  Number 

Probability  density  function 
Observed  wave  height 
Integer  (subscript) 

Maximum  value  of  subscript  1 
Integer  (subscript) 

Maximum  value  of  subscript  j 
Kinematic  (subscript) 

Wavelength 

Total  number  of  observations 

Number  of  observations  in  a specified  cell 

Probability 

Probability  of  slamming 
Relative  sea  state 

Amplitude  of  relative  motion  in  regular  waves 

Significant  single  amplitude  of  relative  motion  in 
random  waves 

Significant  single  amplitude  of  relative  velocity 
in  random  waves 
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Wave  spectral  density 

Significant  single  amplitude  of  absolute  vertical 
acceleration  in  random  waves 

Ship  draft 

Modal  wave  period 

Observed  wave  period 

Ship  speed 

Ship  longitudinal  axis 
A particular  value  (subscript) 

Amplitude  of  regular  wave 
Bow  wave  profile  elevation 
Significant  height  of  random  waves 
Wavelength  to  ship  length  ratio 
Variance  of  random  wave  elevation 
Wave  frequency 
Wave  encounter  frequency 
Modal  wave  frequency 


ABSTRACT 


Methodologies  for  empirical  assessment  of  the  non- 
kinematic  components  of  ship-to-wave  relative  motion,  change 
of  level,  and  bow  wave  profile  are  presented.  Computation 
of  wave  contours  from  observed  wave  data  is  described.  Lin- 
ear superposition  computations  in  the  wavelength/ship  length 
domain  are  outlined. 


ADMINISTRATIVE  INFORMATION 

The  work  described  herein  was  funded  by  the  Conventional  Ship  Sea- 
keeping Research  and  Development  Program  under  Project  Number  62543  and 
Block  Number  SF  43  421  202.  The  work  was  performed  at  the  David  W.  Taylor 
Naval  Ship  Research  and  Development  Center.  There  it  was  identified  by 
Segment  20  of  Work  Unit  Number  1-1504-100. 

INTRODUCTION 

This  document  provides  supporting  details  for  the  1979  Society  of 
Naval  Architects  and  Marine  Engineers  STAR  Symposium  paper  "Minimum  Free- 
board Requirements  for  Dry  Foredecks:  A Design  Procedure,"  Reference  1.* 
The  topics  Involved  are: 

1.  Methodology  for  empirical  estimation  of  nonklnematic  relative 
motion,  change  of  level,  and  bow  wave  profile; 

2.  Derivation  of  wave  environment  characteristics  from  observed 
wave  data;  and 

3.  Linear  superposition  and  related  computations. 

Each  of  these  topics  is  treated  hereinafter. 

NONKINEMATIC  RELATIVE  MOTION,  CHANGE  OF  LEVEL  AND  BOW  WAVE  PROFILE 

Here  the  empirical  formulations  used  to  account  for  the  nonklnematic 
components  of  relative  motion,  for  change  of  level,  and  for  bow  wave  pro- 
file are  described. 


A complete  listing  of  references  is  given  on  page  11. 


RELATIVE  MOTION 

Data  reported  in  References  2 and  3 and  data  from  other  experiments 
which  have  not  been  reported  were  used.  Various  plottings  of  the  ratio 
of  measured  to  kinematic  relative  motion  transfer  function  were  tried. 
(Since  undistorted  wave  dimensions  were  used  on  both  sides  of  the  compari- 
son, this  is  equivalent  to  a direct  comparison  of  measured  to  kinematic 

relative  motion.)  Strong  trends  were  found  with  L /L  and  with  Ship 

w PP 

Station,  but  the  effect  of  Fn  was  hard  to  discern.  In  this  context,  it 

should  be  pointed  out  that  most  of  the  available  data  was  for  Fn  in  the 

0.20  to  0.40  range.  Some  data  applicable  to  Fn  = 0.05  was,  however, 

located.  It  indicated  that  the  ratio  of  measured  to  kinematic  relative 

motion  transfer  functions  was  close  to  unity. 

The  available  data  also  indicated  that  the  ratio  tended  to  unity  as 

L / L increased  to  the  order  of  4.0  and  as  Ship  Station  6 was  approached, 

w pp 

Aft  of  Ship  Station  6,  the  ratio  again  diverged  from  unity,  but  this 

location  was  considered  to  be  a reasonable  stopping  point  for  above-water 

bow  considerations.  No  limiting  trends  in  short  waves  could  be  determined. 

Fortunately,  the  relative  motions  associated  with  short  waves  are  small 

(the  transfer  function  converges  to  unity  in  the  kinematic  case) , so  it 

was  felt  that  the  short-wave  end  of  the  ratio  could  be  handled  in  the 

course  of  spectral  closure  considerations. 

The  available  data  were  ultimately  faired  in  coordinates  of  (r  It  ) / 

a a c 

(r  It,  ),  versus  L /L  with  due  attention  to  the  various  limiting  trends 
a a K w pp 

involved.  A sample  plot  and  fairing  is  shown  in  Figure  1.  Similar  plots 
and  fairings  were  made  for  all  ship  locations  at  which  data  were  available. 
Cross  fairings  over  Ship  Station  were  also  performed.  The  family  of  curves 
presented  In  Figure  2 resulted. 

CHANGE  OF  LEVEL 

Change  of  level  was  evaluated  on  the  basis  of  calm  water  measurements 
for  a number  of  ships.  It  was  found  that  these  measurements  collapsed 
rather  well  at  fixed  Fn  if  ships  with  and  without  bow  domes  were  separated. 
Hence,  it  was  decided  that  a viable  approximation  could  be  obtained  by 
Froude  scaling  from  a similar  hull.  It  is  obviously  important  that  the 
presence  or  absence  of  a bow  dome  be  taken  as  a criterion  for  similarity. 
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BOW  WAVE  PROFILE 

Calm  water  measurements  for  a number  of  ships  were  again  employed. 

As  a rule,  these  measurements  were  available  only  for  design  speed.  In 
those  cases  for  which  such  measurements  were  available  at  two  or  more 
speeds,  it  was  evident  that  the  crest  of  the  bow  wave  moved  forward  as 
well  as  decreasing  in  amplitude  with  decreasing  speed.  It  was  also  evi- 
dent that  the  presence  of  a bow  dome  tended  to  "sharpen"  and  increase  the 
height  of  the  bow  wave's  crest. 

The  latter  factor  simply  implied  that  it  would  again  be  necessary  to 
segregate  ships  with  and  without  bow  domes.  The  shift  in  crest  location 
with  speed,  though,  introduced  a considerable  difficulty.  It  did  not 
appear  that  this  phenomenon  could  be  accounted  for  by  any  simple  scaling 
or  nondimensionalization  scheme. 

To  accommodate  the  shift  in  crest  location  (and  simultaneously  the 
change  in  crest  height  for  a given  ship)  with  speed,  simplified  versions 
of  the  formulas  presented  in  Reference  A were  used.  Specifically,  it  was 
found  that  a point  on  the  bow  wave  profile  at  ship  speed  V^,  say 
where  is  the  wave  profile  height  measured  positive  upward  from  the 
waterline  and  Xq  is  the  ship  location  measured  positive  aft  of  the  forward 
perpendicular,  could  be  transformed  to  approximate  the  bow  wave  profile 
point  at  ship  speed  V using 

S>(x)  = (^)5/4  WV  (1) 


X = (f->  /2  x (2) 

0 

By  applying  equations  (1)  and  (2)  for  a sufficient  number  of  xq's»  then, 
a bow  wave  profile  at  speed  V can  be  constructed. 

It  was  found  that  preceding  the  transformation  defined  by  (1)  and  (2) 
with  Froude  scaling  to  a common  ship  length  caused  the  available  bow  wave 
profile  to  collapse  rather  well  for  subsets  with  and  without  bow  domes. 
Ultimately  then,  the  procedure  for  bow  wave  profile  estimation  involved 
both  Froude  scaling  and  transforming  the  data  for  a similar  hull.  Figure 
3 illustrates  the  procedure  for  a typical  case. 
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DERIVATION  OF  WAVE  ENVIRONMENT  CHARACTERISTICS  FROM  HOGBEN  AND  LUMB  DATA 
The  Hogben  and  Lumb  atlas'*  presents,  for  certain  ocean  areas  and 
seasons,  the  number  of  observations  reported  for  coded  combinations  of 
observed  wave  height,  (bw)0BS»  and  observed  wave  period,  (TW)0BS*  *t  is 
assumed  that  areas  and  seasons  can  be  combined  by  direct  summation  of  the 
relevant  tabulations.  Let  the  result  of  such  a summation  be  represented 

by  ^w^OBS’  ^Tw^OBS  and  "ij  where  1 * 1»  2 1 and  j * 1.  2 

J.  To  each  combination  of  the  subscripts  i and  j there  corresponds  an 

area,  a.,  in  (h  (T  ) space.  The  joint  probability  density 

lj  W UDo  W Udo 

function  of  (hw)0Bg  and  (T^,  say  f{  [ (h^^lj,  [ (V0BS] j > can  thus  be 
empirically  approximated  for  each  (i,j)  combination  by 

f^(hw^OBS^i,^Tw)OBS^j}  “ N a^ 
where  N is  the  total  number  of  observations,  i.e., 

N - £ £ n (4) 

j»l  1*1 


Now,  let  f qC  t ^v^OBS^i’ ^ ^Tw^OBS^ j ^ be  an  arbitrarily  specified  value 
of  f{[(hw)OBSJi‘[(Tw)OBS1jK  Compute 

P0  = ^ ^ aij[f{^(hw)OBS^i,^Tw)OBS^j^] 

(5) 

f0^  ^w^OBsV  ^TwJOBS^  j * 


setting  td^Qgglj}  “ fo{  ^hw^OBS^i’^Tw^OBS^j^  equal  to 

zero  for  all  combinations  of  i and  j such  that  f^{ I ^bw^OBS^i* ^Tw^OBS^j ^ — 

^\p0BS^i’  ^Tw^OBS^j  Tben  the  boundarles  of  tbe  ^hw^OBS*  ^Tw^OBS^ 

cells  for  which  f { [ (hw)0BSli» t (^OBS^ j * “ f0* ^hw^OBS^i’ ^(Tw^OBS^ j * is 

nonzero  roughly  define  a contour  along  which  f{[(h  is 

w Uoa  l w Uoo  j 

constant  at  the  fQ{ [(bw^OBS^i’ ^Tw^OBS^ j * level‘  Further,  PQ  approximates 
the  probability  that  a given  t (bw^oBS* ^Tw^OBS^  observation  will  fall 
within  this  contour. 

One  could,  of  course,  define  any  number  of  surfaces  in  t(bw)0Bg» 
(Tw)qbs^  that  would  in  some  manner  bound  the  fraction  Pq  of  all 
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observations.  Establishing  the  surface  at  a constant  fn[(h  )Anc>(T  )__e] 
is  a rather  arbitrary  expedient.  The  author  finds  this  definition  to  be 
intuitively  appealing.  It  would,  however,  be  of  interest  to  explore 
alternative  definitions. 

We  can,  in  a crude  sense,  identify  the  largest  value  of  f{[(h 

W U do  X 

l(VoBS]j}  for  3 given  value  of  say  max  f{I(hw)OBS1ilI(Tw>OBS1j-m} 

with  the  most  probable  value  of  (h  given  the  specified  (T 

W OBo  W OBf> 

Similarly,  the  most  probable  value  of  (T  ) for  a given  (h  can  be 

identified  by  max  f {[  ('^w)0Bg]  j ) t (hvPoBS^i=m^‘  These  "most  probable" 
values,  together  with  the  constant  f _ [ (h  ),vnc«(T  )r.T>el  contours  and 
associated  probabilities  constitute  the  basic  results  needed  here.  By 
manual  fairing  and  then  a transformation  to  significant  wave  height  versus 
modal  wave  period  coordinates  using  the  Nordenstrom  calibrations^,  wave 
environment  characterizations  such  as  that  presented  in  Figure  3 of 
Reference  1. 


LINEAR  SUPERPOSITION  AND  RELATED  COMPUTATIONS 
Here  several  formulas  presented  in  the  course  of  the  discussion  of 
ship  response  statistic  computations  given  in  Reference  1 are  derived. 

WAVE  SPECTRA 

The  Bretschneider7  wave  spectral  family  can  be  written 

Sc((u)  - Au”  5 e~B<J  (6) 

for  to  from  0 to  infinity.  The  area  under  this  spectrum  is  the  variance  of 
wave  elevation,  i.e., 

00 

\ = / V“)d“ = « (7) 

o 

Further,  the  spectrum  has  a unique  mode  where 

X~fS  (u)]  = 4Bw  - 5 = 0 (8) 

aw  t 
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- . ■: . _ 


. 


Letting  oj  represent  the  mode,  (8)  shows  it  to  be 


“0  (T> 


(9) 


The  wave  period  corresponding  to  toQ  is  taken  to  be  the  modal  period, 
<V  o ’ Thus , 


<Vo  = 2it/<T> 


(10) 


Now  consider  transformation  to  the  L -domain  in  accord  with  the 

w 

gravity  wave  relationship 


“ ■ <£*>* 
w 


(11) 


Under  this  transformation,  (6)  becomes 


AL 


BL  2 
w 


W 2(2*s)2 


W (2irg)2 


(12) 


Integrating  over  S^(LW)  from  0 to  infinity  reproduces  the  right-hand  side 
of  (7).  Differentiating  S (L  ) with  respect  to  L and  equating  the  result 

S'*  W 

to  zero  yields  the  modal  wavelength,  (L  ) , as 

w 0 


2irg 

/2B 


(13) 


The  unit  (£  )i/  condition  can  be  imposed  via  (7): 
w / 3 


1 


Substituting  A and  B from  (13)  and  (14)  into  (12)  yields 


(14) 


W 


w 


16[(Lw)0]2  e 
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l i is; so.:;' - 


(15) 


Nondimensionalizing  the  range  of  (15)  by 


A 


L / L 


v pp 


(16) 


yields 


S?(A) 


XL  2 

_2£ e 

“iW 


Now,  by  defining  the  relative  sea  state  parameter,  R,  to  be 


R 


<Vo 


pp 


equation  1 7 can  be  written 


(17) 


(18) 


S^(X) 


- ~(-)2 
2vr' 


which  is  the  form  given  in  Reference  1. 


MODAL  PERIOD  TO  LENGTH  RELATIONSHIP 

It  follows  from  (10)  and  (11)  that 


<Vo 


as  given  in  Reference  1. 


(19) 


(20) 


NUMERICAL  INTEGRATION 

The  Integrals  included  in  Equations  (7)  through  (1)  of  Reference  1 
must  be  numerically  approximated.  For  those  points  defirtu  by  the  strip 
theory  computations,  this  is  accomplished  by  a straightforward  Lagrangian 
integration  procedure.  The  difficult  part  of  the  problem  is  approximating 
the  "tails”  of  the  response  spectra  which  lie  beyond  the  computed  points, 
i.e.,  obtaining  spectral  closure.  In  the  cases  of  relative  motion  and 
velocity,  which  are  of  major  concern  here,  the  problem  is  emphasized  by  the 


limiting  behavior  of  the  associated  response  amplitude  operators.  As  u>e 

becomes  large,  [ (r  /?  ),  ]2  approaches  unity  and  [ (r  /?  ),  ]2  (=  w 2[(r  / 

^ a a tc  a a k e a 

C ).]Z)  becomes  very  large. 

SL  k 

Working  in  the  L /L  domain  minimizes  the  problems  associated  with 
w pp  F 

the  high  frequency  behavior  of  relative  motion  and  velocity.  In  the  u)e 

domain  this  behavior  theoretically  continues  to  infinity.  In  the  L^/L^ 

domain,  it  is  confined  to  a small  range  near  zero. 

With  L /L  identified  by  A,  let  min(A)  be  the  smallest  value  of  A 
w pp 

for  which  strip  theory  computations  are  performed.  Then  [(r  /c  and 

[ (ra/ c 3 2 are  extrapolated  from  their  computed  values  at  min(A)  to  unity 
at  A = 0.  The  particular  values  of  A for  which  extrapolated  values  are 
computed  are  0.01,  min(A)/4,  min(A)/2  and  3[min(A)]/4.  The  extrapolated 
points  are  included  in  the  Lagrangian  integration. 

Low  A closure  for  the  acceleration  integral  is  handled  by  extra- 
polation at  the  response  spectrum  level.  It  is  assumed  that  the  accelera- 
tion spectrum  varies  in  a linear  manner  from  its  computed  value  at  min(A) 
to  zero  at  A = 0.  The  resultant  triangular  area  is  included  in  the 
spectral  area. 

High  A closure  is  handled  the  same  for  all  responses.  This  is  done 
at  the  spectral  level.  A least  squares  line  is  fitted  to  the  response 
spectrum  ordinates  at  the  three  largest  values  of  A for  which  PHM  computa- 
tions are  performed.  The  area  bounded  by  this  line,  by  the  A axis,  and 
by  the  spectral  ordinate  at  the  largest  value  A is  then  included  in  the 
spectral  area.  If  the  slope  of  the  fitted  line  is  positive,  indicating 
that  closure  is  not  attained,  a flag  is  set.  This  has  not  occurred  in  any 
computations  performed  to  date.  In  fact,  the  contribution  of  the  high  A 
closure  area  to  response  per  unit  wave  height  was  generally  insignificant. 

LIMITING  WAVE  HEIGHT  FOR  SLAMMING 

The  probability  of  bottom  slamming  from  Reference  8 can  be  written 

■ 2[(F?")2  + (f-1-)2! 

P - e (21) 

s 

In  terms  of  response  per  unit  wave  height,  (21)  can  be  rewritten  in  the 
form 
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r 


i 


in  that  the 
specifying  a 
Uw)l/3.  say 

(24) 

Equation  (12) 


(This  equation  differs  from  Equation  (13)  of  Reference  1 only 

notation  of  the  latter  specializes  it  to  Ship  Station  3.)  By 

numerical  value  for  P , say  lim(P  ),  a corresponding  value  of 

s s 

lim[(£  )i,  ]CTA„,  can  be  derived  from  (22).  It  is 
W / 3 oLAM 

l±mr (Cw) !/  ]SLAM  = *ln[iim(p  )]}* 

s 

which,  again  except  for  specialized  notation,  is  identical  to 
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Figure  1 - Relative  Motion  Transfer  Function  Ratios 
in  Kay  of  Ship  Station  0 
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Figure  2 - Composite  to  Kinematic  Relative  Motion 
Transfer  Function  Ratios 
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Figure  3 - Bow  Wave  Profile  Scaling  and  Transformation 
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THEY  CARRY  A DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 
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MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A CASE-BY  CASE 
BASIS. 


